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ABSTRACT 

We develop a statistical method for using multicolor photometry to 
determine distances using Cepheid variables including the effects of temperature, 
extinction, and metallicity and apply it to UBVRIJHK photometry of 694 
Cepheids in 17 galaxies. We derive homogeneous distance, extinction and 
uncertainty estimates for four models, starting from the standard extragalactic 
method and then adding the physical effects of temperature distributions, 
extinction distributions, requiring positive definite extinctions, and metallicity. 
While we find general agreement with published distances when we make similar 
systematic assumptions, there is a clear problem in the standard distances 
because they require Cepheids with negative extinctions, particularly in low 
metallicity galaxies, unless the mean LMC extinction exceeds E(B — V) ^ 0.25. 
The problem can be explained by the physically expected metallicity dependence 
of the Cepheid distance scale, where metal-poor Cepheids are hotter and fainter 
than metal-rich Cepheids. For V and I we found that the mean magnitude 
change is —0.14 ± 0.14 mag/dex and the mean color change is 0.13 ± 0.04 
mag/dex, with the change in color dominating the change in distance. The effect 
on Type la supernova estimates of the Hubble constant is dramatic because 
most were found in the metal poor galaxies with the bluest Cepheids. The Type 
la Multi-color Light Curve Shape (MLCS) method estimate for H formally 
rises from 69 ± 8 km s _1 Mpc -1 to 80 ± 6 km s _1 Mpc -1 with the metallicity 
correction. 

Subject headings: Cepheids - distance scale - galaxies: distances and redshifts - 
Hubble Constant 
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1. Introduction 

Cepheid variables are fundamental to most extragalactic distance estimates, 
determinations of the Hubble constant, and models for the structure of the Galaxy because 
the Cepheid period-luminosity (PL) relations are accepted as one of the most accurate 
primary distance indicators. In the last few years, the number of extragalactic Cepheids 
has exploded due to the two large extragalactic surveys using the Hubble Space Telescope 
(the Extragalactic Distance Scale Key Project (Freedman et al. 1994, Ferrarese et al. 1996, 
Graham et al. 1996, Kelson et al. 1996, Silbermann et al. 1996) and the Type la Supernova 
Calibration Project (Saha et al. 1994, 1995, 1996ab, 1997)), and the microlensing surveys of 
the Magellanic Clouds (e.g. Beaulieu et al. 1996, Welch et al. 1997). The goal of the HST 
projects is to determine the Hubble constant with 10% (0.2 mag) accuracy, which requires 
uncertainties in the Cepheid distance estimates, including all systematic uncertainties, that 
are still smaller. In order to reach these goals, the observational projects (see the review 
by Freedman 1997) have expanded the range of filters used to study Cepheids, particularly 
into the infrared, introduced systematic corrections for extinction, and tried to find good 
empirical tests for the effects of metallicity on Cepheid distance estimates. 

The physical basis for Cepheids as distance estimators rests on two foundations (see the 
reviews by Feast & Walker (1987), Madore & Freedman (1991), and Tanvir (1996)). The 
first foundation is stellar structure, which closely correlates the mass, radius, luminosity, 
temperature, and oscillation period of a star, so that in any localized region of the H-R 
diagram there is period-luminosity-color ("PLC") relation stating that the luminosity can 
be determined from the color (temperature) and the period (mass and radius) of oscillation. 
The second foundation is that the physics of the oscillations limits the Cepheids to a narrow 
range of luminosity and temperature called the instability strip. As a result, any projection 
of the three-dimensional PLC space onto a two-dimensional subspace produces a tightly 
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correlated relation. In particular, the period-luminosity (PL) correlations are defined by 
projecting over the temperature/color distribution. Both the PLC relations and the location 
of the instability strip are expected to be functions of composition (e.g. Stothers 1988, Stift 
1990, Chiosi, Wood & Capitanio 1993), although there is great debate about the magnitude 
of the dependence and its measurement (e.g. Caldwell & Coulson (1986, 1987), Freedman 
& Madore (1990) Gould (1994), Stift (1995), Sasselov et al. (1996)). The sense, however, is 
that metal rich Cepheids are cooler and brighter than metal poor Cepheids at fixed period. 

The goal of any analysis of Cepheid data is to accurately determine the distance to 
the Cepheid and its uncertainty, after compensating for the effects of period, temperature, 
composition, and extinction. Existing Cepheid analysis methods are divided into 
extragalactic and Galactic approaches. For extragalactic systems the analysis must 
determine the common distance of an ensemble of Cepheids, usually based on two-color, 
low-accuracy photometry, with poor phase coverage. The standard approach (see Madore & 
Freedman 1991) uses only the correlation between luminosity and period (the PL relations) 
to estimate the distance modulus and the mean extinction. Gould (1994) pointed out 
that the method is statistically inefficient because it ignores the strong correlations in the 
residuals (i.e. the PLC relations). Galactic analyses (e.g. Caldwell & Coulson 1986, 1987, 
Caldwell & Laney 1991, Laney & Stobie 1993, 1994) must estimate the distances and 
extinctions of individual Cepheids using three-color photometry and PLC relations because 
the Cepheids no longer lie at a common distance. The standard extragalactic approach 
explicitly ignores the PLC relations, and the standard Galactic approach (in some senses) 
ignores the instability strip. 

In the following sections we develop a self-consistent physical and mathematical 
analysis of multi-color Cepheid mean magnitudes as the first in a series papers reanalyzing 
the Cepheid distance scale. Our approach differs from the standard approaches in three 
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major ways. The first difference is that we analyze all the data simultaneously rather than 
one galaxy at a time. The use of Cepheids as distance indicators is predicated on the 
homogeneity of their physical properties, so either all the data can be self-consistently and 
simultaneously analyzed or we must reject the Cepheids as distance indicators. Moreover, 
the distance estimates (and other physical parameters) are very highly correlated. Any 
change in the distance or extinction estimate for one Cepheid galaxy requires a simultaneous 
change in the values for all Cepheid galaxies. The second difference is that we avoid the false 
dichotomy between the Galactic and extragalactic analysis procedures, and demonstrate 
how to reconcile the two approaches. The resulting statistical method has greater statistical 
efficiency than traditional methods, can model all measured colors simultaneously, and 
allows for better control, treatment, and understanding of systematic problems such as 
individual Cepheid extinctions and the effects of metallicity. The resulting scheme, although 
based on a physical model, closely resembles the empirical treatments of Gould (1994) and 
Sasselov et al. (1996). The third difference is that we explore the important systematic 
errors that can affect the Cepheid distance estimates, particularly the physics of extinction, 
temperature and metallicity. While treatments of extragalactic Cepheid distances generally 
recognize the existence of these systematic uncertainties, they are rarely included in the 
distance estimates or their uncertainties. We develop our analysis method in §2, and 
compare it to the existing techniques. In §3 we examine the Cepheid correlations (PL, 
PLC relations etc.) and the resulting distance estimates. We summarize our results, their 
shortcomings, and possible solutions in §4. 

2. Cepheid Distance Estimates 

We assume that the bolometric luminosity and the bolometric correction are linear 
functions of the period p = log P/Pq, the effective temperature t = log T e /T , and the 
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( logarithmic) metallicity Z, so that the (intensity mean) apparent magnitude V k of a 
Cepheid in band k, at distance modulus /i, with extinction E is 

V k = V ok + a k p + (3 k t + lk Z + fid k + ER k (1) 

where V ok is the magnitude zero-point vector, R k is the reddening vector, d k = 1 is the 
distance vector, and logPo — 1-4 and T are reference periods and temperatures. We 
assume that the slopes a k and f3 k are independent of metallicity, although theory predicts 
a weak dependence (see Stothers 1988, Stift 1990, 1995, Chiosi et al. 1993). If we neglect 
the dependence of bolometric corrections on surface gravity, a k = otj. We call eqn. (1) 
a period-luminosity-color (PLC) relation, although standard PLC relations (see Feast & 
Walker 1987) replace the effective temperature by a color. We neglect additional variables 
such as the Cepheid's age, the instability strip crossing number, the helium abundance, and 
non-linearities in the PLC relation (see Caldwell & Coulson 1986), so we assume that the 
relations defined by eqn. (1) have an intrinsic width of o PLCk - 

Formally, if we know the precise values of the vectors a k , f3 k) 7^ and R k , the vectors are 
not degenerate with d k (= distance), and we possess accurate photometry of the Cepheid 
in a sufficiently large number of bands, then the PLC relations can be used to determine 
the distance to a particular Cepheid. Unfortunately, since we must determine a k , j3 k and 
7fc as we proceed from a small number of colors, it is difficult to use the PLC relations to 
determine distances without adding additional constraints. The additional constraints used 
in standard Cepheid analyses (e.g. Madore & Freedman 1991) correspond to adding priors 
on the distances and extinctions. We now develop a general mathematical description for 
fitting Cepheid magnitudes that includes the standard methods as subcases or limits of a 
more general model. 

The PLC relations contain no information about which stars pulsate, and the most 
important prior information is the location and width of the instability strip. The 
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standard extragalactic method uses the instability strip by averaging the distribution over 
temperature and using only the PL relations to determine distances. We parametrize the 
instability strip by a period-color (PC) relation defining the temperature distribution for 
stars at a fixed period using a likelihood function 



L(t\p) oc exp 



(t - 5 Z Z - S P p) 
2a\ 



(2) 



'PC 

including a possible shift in the location of the instability strip with metallicity (e.g. Stift 
1990, Chiosi et al. 1993). We assume that the width of the instability strip is constant, 
although the observed narrowing of the strip at short periods (e.g. Fernie 1990) could be 
modeled by making ape a function of period (e.g. Gould 1994). We chose to parametrize 
the instability strip with a PC relation rather than a luminosity-color relation (i.e. the 
H-R diagram) because it is distance independent and uses the well-defined period as the 
independent variable instead of the luminosity. From these two assumptions we can derive 
all the standard relations used to study Cepheids and their correlations. For example, the 
PL relation in band V k is 

T^wf = {Vk) = v ° k + a ' kp + ii{z) + ^ + {E)Rk (3) 

where a' k = ak + fikSp and 7^ = 7fc + Pk^z, and the dispersion in the relation is 
a PLCk + Pk a PC — Pl a PC smce a PLCk <^ °"pc- From here on we use the deviation of the 
temperature from the expected mean, St — t — SzZ — Spp, which also shifts the period 
and metallicity vectors to the a' k and 7^ appearing in the PL relation (3). We cannot 
independently determine a k , 7/c, Sz and Sp without absolute temperature references, so we 
restrict our solutions to determinations of St, a' k , and 7^. We measure distances, extinctions, 
and metallicities relative to the LMC values of {[i)lmc, (E)lmc, and Zlmc = [0/ H]lmc- 
We use the oxygen abundance for the metallicity variable because it is the only abundance 
available for most of the extragalactic systems. We can also constrain the models using prior 
information on the distance modulus, extinction, and metallicity. In external galaxies we 
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can use the constraint that all Cepheids lie at the same distance, but this assumption begins 
to fail for the LMC and may be a poor assumption for the SMC (e.g. Caldwell & Coulson 
1986, Caldwell & Laney 1991). We allowed the Magellanic Clouds to be tilted relative to 
the line of sight, and to have a finite thickness. Distances and extinctions were fit for each 
Galactic Cepheid, constrained by the non-linear likelihood of fitting the measured radial 
velocity with a flat rotation curve. From the rotation curve model we obtain estimates 
of the solar radius Ro and circular velocity ©o- We model the priors using Gaussians 
with mean values of (fx) for the distance modulus, (E) for the extinction, and (Z) for the 
metallicity, with dispersions of <r M , erg, and az for each galaxy or Cepheid. The relative 
calibration of the standard and the HST V and I magnitudes is uncertain at the level of 
0.05 mag (Hughes et al. 1994), so we include two HST calibration variables constrained 
by a Gaussian prior of width 0.05 magnitudes in the likelihood. The temperature ape and 
PLC c PLCk widths were the same for all galaxies. 

We can divide our analysis into two parts. First, we estimate the parameters of the 
individual Cepheids (/i, 5t, E, and Z) given the current parameters of the global model. 
Second, we optimize the parameters (Vofc, o.' k1 7^,, • • •) of the global model. For a 
particular Cepheid we have mean magnitude measurements V m k with uncertainties a m k in 
each of k — 1 • • • N bands and a known period P. If we define a\ = <7^ fe + <Jp LCk , then 
the log-likelihood for the model to fit the measured mean magnitudes V m k of a particular 
Cepheid is 

- 2 InL = £ V* ~ YH + 41 + < £ - f» a + it ~ + (g - f>? + ln 1 s -.| (4) 

fe=l ^ °PC °E °l °Z 

up to a constant, where the covariance matrix is 

S' 1 = + vldidj + PiPjtrpc + RiRjol + 7-7^1 > ( 5 ) 
and the indices run over the filters included in the calculation. We must use the determinant 
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of the covariance matrix S in the likelihood if we are to simultaneously determine the 
properties of the individual Cepheids and their global statistical relations. 

We can relate our model to standard analyses by breaking the calculation into two 
sections: first, deriving the deviations in the properties of individual Cepheids from the 
mean, and second, determining the PLC relation vectors and the mean properties of the 
Cepheids. If we measure the magnitude residuals A 14 = V m k — (14) relative to the mean PL 
relations (eqn. (3)) for the parent galaxy, define the deviation of the Cepheid parameters 
from the mean parameters by x and the vector- weighted residuals by v, 



( 



x 



St 

E-(E) 

\ Z-(Z) ) 
and define the matrix C by 

°k a PC 
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(7) 



then x = C~ l v and the covariance matrix of the parameter estimates for a fixed PLC 
relation is C^ 1 . With the optimization of these variables, the contribution of the Cepheid 
to the likelihood becomes 

L oc IS]' 1 / 2 exp(--AV T SAV) (8) 
2 

where Ay is the vector of residuals relative to the mean PL relations and S is the covariance 
matrix defined in eqn. (5). 



We obtain the standard extragalactic method (see Madore & Freedman 1991) if the four 
priors have zero width (ape = &e = = cr z = 0) and there is no metallicity dependence 
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il'k = 0)- The covariance matrix S^- = S i j(al li + <J 2 PLC i) is diagonal. Some extragalactic 
Cepheid distances are derived by using reddening free distances or determining individual 
extinctions (e.g. Freedman et al. 1990, 1991, 1992, Tanvir et al. 1995, Saha et al. 1996ab, 
1997), which corresponds to allowing oe — *• oo. If we optimize the widths of the priors, then 
the method matches that of Gould (1994) or Sasselov et al. (1996) in using the covariance 
matrix of the residuals to build a better statistical model of the data. The advantage of our 
formalism is that it derives from a physical model and we gain new physical insights into 
the system from the variables making up the covariance matrix. The disadvantage is that if 
our physical model is incorrect or it is missing important sources of correlated variance in 
the data, it is not as optimal or correct a statistical approach as using a purely empirical 
covariance matrix. If we limit the model to two or three filters and use broad priors we 
recreate the normal Galactic approach (e.g. Caldwell & Coulson 1986, 1987, Laney & 
Stobie 1986, 1993, 1994, Fernie 1990, Fernie et al. 1995) with the color dependence of the 
standard PLC relations appearing indirectly through the temperature variable. 

Three classes of data and parameters enter the problem. The first class consists of 
the period and the composition. Here we know the dependent variable (p or Z), and seek 
to determine the state vectors (a' k and 7^). Given an adequate range for the dependent 
variables in an external galaxy, the state vectors are well-determined and non-degenerate. 
In particular, the standard PL relations (e.g. Madore & Freedman 1991, Laney & Stobie 
1994, Tanvir 1996) determine a' k based on the Cepheids in the LMC. In the second class, 
consisting of the extinction and the distance, we know the state vectors (R k and d k ) and 
would like to determine the dependent variable (E and p). The distance vector is known 
exactly (d k = 1), and the extinction vector R k is known approximately (see §2.2). With 
accurate measurements in a sufficient number of bands, we can determine the distances 
and extinctions for individual Cepheids. The uncertainties will depend on how well we 
can determine the effects of temperature and metallicity, but the formalism will correctly 
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include these uncertainties in the distance and extinction estimates. 

In the third case, the temperature, we know neither the dependent variable t nor the 
state vector (3k- One immediate consequence is the existence of a mathematical degeneracy 
under a rescaling of the temperature by t — > £t. The likelihood is unchanged if we rescale the 
other temperature related variables by (3k — > t^^fik, &p £8p, $z —> £8z and a PC — > £crpc- 
No observable (i.e. magnitude) depends on the rescaling, and it means that we cannot 
set an absolute temperature scale. A more fundamental problem is that values found for 
variables such as the t—/3k pair may not have physical meaning assigned to them in the 
mathematical model. The minimization procedure will simply use them to absorb as much 
variance as possible from the residuals, and if (3k is unconstrained it assumes the value of 
the most important unmodeled principal component of the true covariance matrix. Only 
if the primary source of the variance is the temperature distribution at fixed period will t 
represent the physical temperature. In essence, we are defining the covariance matrix in 
terms of a few principal components defined by the state vectors, and we assert that the (3k 
principal component represents temperature. The degeneracy only affects our interpretation 
of the variables, and the statistical model will still reproduce the correct observational 
PLC/PC/PL relations. For practical purposes we solved the degeneracy problem by using 
a prior for the (3k derived from fitting Cepheid light curves, which we discuss in paper II 
(Kochanek 1997b). The difficulty in estimating (3k (or its equivalent slope in standard PLC 
relations) and separating temperature from extinction has lead the extragalactic Cepheid 
community to avoid modeling the temperature distribution (e.g. see the critique of PLC 
relations in Madore & Freedman 1991). 
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2.1. Data 

We simultaneously analyzed the Cepheids of 17 galaxies (see Table 1), including only 
Cepheids with periods between 7 days and 80 days. We used all available Johnson UBV, 
Kron-Cousins RI, and Glass-Carter JHK mean magnitudes for the Cepheids, because the 
systems with many color photometry (the Galaxy LMC, SMC, M 31, M 33, and NGC 300) 
offer the best hope of separating the effects of distance, temperature, extinction, and 
metallicity. With two-color photometry there is no simple means of separating the effects 
of temperature and extinction, while with eight-color photometry it should be possible. We 
assigned magnitude uncertainties of 0.05 mag for the Galactic, LMC and SMC Cepheids, 
and 0.10 mag for the HST Cepheids and the poorly sampled ground-based data on M 31, 
M 33, and NGC 300. Estimates of the distances and extinctions were little affected by 
changes in the estimated measurement errors. Although we spot-checked many of the 
inferred magnitudes and periods for the Cepheids, we used the published periods and 
magnitudes in our analysis. 

Galaxy: We used the Cepheids with radial velocities in Pont et al. (1994) and periods 
longer than 7 days. Intensity mean magnitudes were checked from the data 
compilations by Welch (1996) and Berdnikov (1987, 1996), and matched existing B 
and V tabulations (e.g. Fernie et al. 1995). We dropped the Cepheids AA Ser, XZ 
Car, and SU Cru from the Pont et al. (1994) sample as outliers. 

LMC: We used the UBVI data from Caldwell (1996), the JHK data from Laney k Stobie 
(1986, 1993, 1994) which includes earlier data by Welch et al. (1987), and the R 
data from Madore (1985). The Caldwell (1996) mean magnitudes include earlier data 
whose sources are reviewed in Madore (1985). We rejected HV2301, HV2378 and 
HV2749 from our fits as outliers in the likelihood. Caldwell & Laney (1991) and 
earlier authors also report that these three Cepheids have abnormal properties. 
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SMC: We used the BVRI data from Caldwell (1996) and the JHK data from Laney & 
Stobie (1993, 1994). We rejected HV854, HV1369, HV1438, HV1482, HV1484, and 
H VI 695 as outliers in the likelihood. Caldwell & Laney (1991) and earlier authors 
report HV1369, HV1484, HV1636, and HV1641 as outliers. We find nothing peculiar 
about HV1636, and HV1641 was not included in the sample from Caldwell (1996). 

M 33: We use the BVRI data from Freedman, Wilson & Madore (1991) excluding 21979, 
23764, and Bl. The photometry for V12 in Table 5 of Freedman et al. (1991) is 
shifted to the left by one column. None of the remaining Cepheids stood out in the 
likelihood distribution. 

M 31: We use the BVRI data from Freedman & Madore (1990, Freedman 1996). We 
rejected the 18.5 day period Cepheid in Baade's Field III. 

NGC 300: We use the BVRI data from Freedman et al. (1992). We rejected V21 as an 
outlier. Freedman et al. (1992) comment that the light curve of V21 is not well 
determined. 

M 81: We use the VI data from Freedman et al. (1994). We left the magnitude calibrations 
unchanged. We rejected C8 as an outlier. 

M 100: We use the VI data from Farrarese et al. (1996). There is a typographical error in 
the distance modulus given in the abstract of Farrarese et al. (1996, erratum 1997), 
where the correct value is /i = 31.04 ± 0.17. We added the 0.05 mag "long exposure 
correction" to the tabulated Cepheid magnitudes (Hughes et al. 1994). We rejected 
the two shortest period Cepheids, C68 and C70 as outliers. 

M 101: We use the VI data from Kelson et al. (1996) including all the I photometry (the 
"weak photometry restriction"). The light curve data tables of Kelson et al. (1996) 
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have the time and magnitude columns in different orders. We rejected C28 as an 
outlier. 

IC 4f82: We use the VI data from Saha et al. (1994). We rejected C3-V9 as an outlier. 

NGC 5253: We use the VI data from Saha et al. (1995). We rejected C3-V1 and C3-V2 as 
outliers. 

NGC 4536: We use the VI data from Saha et al. (1996a). Following the authors we 
included only Cepheids with periods longer than 20 days and quality class 4-6 in 
the standard analysis. We added the 0.05 mag "long exposure correction" to the 
tabulated Cepheid magnitudes. We rejected C3-V12, C3-V16, C3-V24, and C3-V31 
as outliers. 

NGC 925: We use the VI data from Silberman et al. (1996). We rejected 3-11, 3-29, 
4-68, and 4-71 as outliers. 

M 96: We use the VI data from Tanvir et al. (1995). 

NGC 3351: We use the VI data from Graham et al. (1996) but restricted our sample to 
periods longer than 10 days. We rejected C7, C17, C23, and C46 as outliers. 

NGC 4496A: We use the VI data from Saha et al. (1996b). Following the authors we 
included only Cepheids with periods longer than 17 days and quality class 4-6 in 
the standard analysis. We added the 0.05 mag "long exposure correction" to the 
tabulated Cepheid magnitudes. 

NGC 4639: We use the VI data from Saha et al. (1997). Following the authors we 
included only Cepheids with periods between 20 and 63 days and quality class 4-6 
in the standard analysis. We added the 0.05 mag "long exposure correction" to the 
tabulated Cepheid magnitudes. 
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2.2. Extinction 

We require an extinction vector R k = A k /E(B — V) defined for all eight niters. Our 
standard vector is R ok = (5.05,4.31,3.30,2.73,2.07,0.95,0.64,0.39) for the UBVRIJHK 
filters based on the Cardelli, Clayton & Mathis (1989) model for the extinction curve. We 
fixed Ry = 3.3 to match the Key Project, but as emphasized by Cardelli et al. (1989) 
the difference between Ry = 3.1 and Ry = 3.3 is mainly in the absolute normalization of 
the R k vector and the estimated E(B — V). The shape of the extinction vector changes 
very little. The Type la project uses A v /Aj = 1.7 instead of A v / ' Ai = 1.6. The earlier 
M 31 (Freedman & Madore 1990) and M 33 (Freedman et al. 1991) studies used R v = 3.1 
combined with the Cardelli et al. (1989) extinction curve, while Ry = 3.3 was used for 
NGC 300 (Freedman et al. 1992). Laney & Stobie (1993) estimated that the JHK values 
were (0.82,0.49,0.30). We do not include the temperature variations in the extinction 
coefficient used by many of the Galactic Cepheid models (e.g. Caldwell & Coulson 1986, 
Laney & Stobie 1986, 1993, 1994, Feast 1987, Fernie 1990). 

Cepheid distances are sensitive to the assumed structure of the extinction vector 
Rk, but most existing treatments of Cepheid distances treat the extinction as a known, 
understood property of the ISM even while different analyses select different models. 
Exceptions are Laney & Stobie (1993, 1994) and Sasselov et al. (1996), who systematically 
varied the R k or tried to determine them from the data, and Freedman et al. (1990, 1991) 
who examined the effects of using standard Ry = 3.1 or 3.3 or extremal models on the 
distance estimates. We will use a fixed extinction model for all the Cepheids, but we 
allow the coefficients to adjust themselves to best fit the data, and the uncertainties in the 
extinction vector will be included in the distance uncertainties. We constrain the extinction 
vector R k to fit the Cardelli et al. (1989) extinction vector R ok by adding a Gaussian prior 
to the likelihood oc exp(— (Rok — Rk) 2 /2<t|j) with or = 0.1 and R y = Ry. The uncertainty 
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roughly matches the uncertainties in the Cardelli et al. (1989) model and the differences 
between various extinction models used for Cepheid distances. 

The extinction is not a free variable because it must be positive definite and larger 
than the estimated foreground extinction (taken from Burstein & Heiles (1984), see Table 
1). Freedman et al. (1990, 1992) were the first to face this problem when they found a 
negative estimated extinction for the Cepheids in M 31 (Baade IV field) and NGC 300, 
although Sasselov et al. (1996) were the first to use the positivity of the extinction as a 
constraint when estimating the metallicity dependence of the Cepheid distances from the 
EROS sample (Beaulieu et al. 1995) of LMC and SMC Cepheids. Most treatments have 
either advocated raising the mean LMC extinction (Freedman et al. 1992, Bohm-Vitense 
1997), ignored the problem (e.g. Saha et al. 1997), or used a statistically incorrect solution 
in setting E — without changing the distance. When we search for solutions with physical 
extinctions, we do so by adding a term to the logarithm of the likelihood (eqn. (4)) of 
the form (E ki — Ef k ) 2 /a% if the extinction E ki of Cepheid i in galaxy k is less than the 
estimated foreground Ef k , with the scale of <je = 0.045 set to include the uncertainty in the 
foreground extinction and the mean LMC extinction. Our choice for the penalty function 
has the advantage of simplicity, but it rises too abruptly to be an ideal model for global 
uncertainties in the extinction scale. The structure of the penalty determines the weight 
attached to Cepheids whose magnitude uncertainties produce negative extinction estimates 
even if the true extinction is positive. We experimented with several more complicated and 
mathematically graceful forms for the penalty, but changes in implementation had little 
effect on the physical results. 
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2.3. Metallicity 

There is no debate about the qualitative effects of metallicity on the Cepheids: 
high metallicity Cepheids are brighter and cooler than low metallicity Cepheids at fixed 
period. The debate centers only on the magnitude of the effect on distance estimates. 
The flux is reduced in the blue due to line blanketing, and increased in the red and 
infrared by backwarming. Theoretical estimates suggest the effect is moderate, about 
(i'b + 7v)/2 — 0.14 mag/dex and (7^ — 7^) = 0.16 mag/dex for the mean B and V band 
luminosity and color (Stothers 1988). Chiosi et al. (1993) find a weaker effect of about 
(l'v + 7/)/2 = —0.05 to 0.13 mag/dex in the mean magnitude, and (7(7 — 7}) = 0.05 to 
0.09 mag/dex in the V — I color. An obvious signature of a metallicity dependence in 
the observational data should be a correlation of Cepheid extinction estimates with the 
metallicity of the parent galaxy. 

Recent attempts to measure the metallicity effects either examined M 31 or compared 
the LMC and SMC. Freedman & Madore (1990) examined the Cepheids in M 31 where the 
work of Blair, Kirshner & Chevalier (1982) implied a steep metallicity gradient. They found 
a distance change with metallicity of 5fi ~ (0.32 ± 0.21) mag/dex where the numerical 
estimate of their gradient is due to Gould (1994). The weak significance of the Freedman & 
Madore (1990) estimate is used by the Key Project as the basis for neglecting the effects of 
metallicity in distance estimates pending an analysis of metal rich and metal poor Cepheids 
in M 101. Gould (1994) also reanalyzed the M 31 data including the empirical covariance 
matrix of the residuals to find a larger correction of (0.88 ± 0.16) mag/dex, although 
the numerical value depended on the colors used. The low metallicity Baade IV field 
required a large negative extinction, suggesting that the metallicity requires a color term, 
as emphasized by Stift (1995) and Sasselov et al. (1996), unless the mean LMC extinction 
is significantly underestimated (Freedman et al. 1992, Bohm-Vitense 1997). Sasselov et 
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al. (1996) analyzed the EROS sample of fundamental and overtone Cepheids in the LMC 
and SMC (Beaulieu et al. 1995) using a method based on Gould (1994) but including 
both distance and color terms for the metallicity. Unlike the M 31 studies, where the 
estimate depends on Cepheids of differing metallicity at a common distance, the Sasselov 
et al. (1996) value relies on positivity of extinction and estimates of the foreground and 
internal extinctions of the LMC and SMC to estimate the effect. They find a correction 
of 0.4+o;2 mag/dex in the mean magnitude and (0.20 ± 0.02) mag/dex in the V-I color. 
The color dependence is similar to the Caldwell & Coulson (1986, 1987) and Gieren et 
al. (1993) estimates for the Galactic Cepheids, of (B — V) oc (0.29 ± 0.05) mag/dex and 
(V-I)oa (0.20 ± 0.05) mag/dex. 

We use the logarithmic abundance of oxygen relative to the LMC (see Table 1, 
AZ = [O/H] — [0/ H] LMC ) as our metallicity variable because it is the only measured 
abundance for most of the Cepheid galaxies. Where possible we include the spatial gradients 
of the metallicity from Zaritsky, Kennicutt & Huchra (1994) and assign metallicities to the 
Cepheids based on their positions in the galaxy (M 33, M 31, NGC 300, M 81, M 101, 
NGC 925, NGC 3351, M 100). Where no gradient was available we used a mean metallicity 
for the whole system (LMC, SMC, IC 4182, NGC 5253, M 96, NGC 4536, NGC 4496A). 
NGC 4536, NGC 4496A, and NGC 4695). We assigned the Galactic Cepheids a metallicity 
0.3 higher than that of the LMC, and we decided not to force a radial metallicity gradient 
in our current model. For consistency with Freedman & Madore (1990) and Gould (1994) 
we used the Blair et al. (1982) gradient for M 31 rather than the Zaritsky et al. (1994) 
gradient. Zaritsky et al. (1994) averaged the Blair et al. (1982) and Dennefeld & Knuth 
(1981) data even though Blair et al. (1982) strongly disagreed with Dennefeld & Knuth's 
(1981) results. We also find an anomalously high metallicity estimate for the M 33 Cepheids 
compared to M 31. Both of these potential problems could lead to an underestimate of the 
strength of any metallicity effect. 
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3. Statistical Models of Cepheid Magnitudes and Distances 

In this section we build the full model for the Cepheid mean magnitudes, starting 
from the standard extragalactic method in Model 0. Since the residuals are roughly 
parallel to the extinction vector, Model 1 allows all Cepheids to have individually estimated 
extinctions. In Model 2 we allow the Cepheids to have a distribution in temperature at 
fixed period, and enforce the positivity of the extinction on the models. Finally, we estimate 
the effects of metallicity in Model 3. We used a flat rotation curve model for the Galaxy to 
obtain estimates for the solar radius R and circular velocity O . The LMC and SMC are 
tilted relative to the line-of-sight, have zero width, and their relative distance is determined 
between the (bar) centers. The tilt is optimized as part of the solution and agrees well 
with Caldwell & Laney (1991). To facilitate comparisons with the standard PL relations of 
Madore & Freedman (1991), we assume an LMC distance modulus of 18.5 mag and a mean 
LMC extinction of E(B — V) = 0.10, but we have shifted the period origin to logP = 1-4. 
We only discuss the results for the Johnson UBV, Kron-Cousins RI, and Glass-Carter JHK 
bands. 

3.1. Model 0: The Standard Method 

We start our analysis by considering the Cepheid correlation functions in the absence 
of metallicity (7^ = 0), scatter in the temperature ((3k = 0), and allow only the Galactic 
Cepheids to have individually determined extinctions. Model treats all the extragalactic 
Cepheids using the standard extragalactic method. The parameters for the PLC relations 
and their uncertainties are summarized in Table 2, and the derived distances and mean 
extinctions for the extragalactic samples are summarized in Table 3. Our zero-points and 
period slopes are in general agreement with Freedman & Madore (1991). The zero-points 
are approximately 0.04 mag and 0.02 mag fainter in the critical V and I bands, but 
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the differences lie well within the standard uncertainties. Our PL relation slopes a' k are 
uniformly shallower by 0.05-0.10, although the agreement is again consistent with the 
uncertainties in the individual estimates. The zero-points and slopes agree with recent 
studies by Laney Sz Stobie (1994)[| and Tanvir (1996) .0 The best fit extinction vector differs 
from the nominal Cardelli et al. (1989) vector by 0.15 ± 0.08, -0.01 ± 0.06, -0.07 ± 0.05, 
-0.05 ± 0.05, +0.10 ± 0.06, +0.02 ± 0.06, and 0.08 ± 0.06 for U, B, R, I, J, H, and K 
respectively. The changes in the coefficients are not a simple change in the Ry value of 
the Cardelli et al. (1989) model, and the final uncertainties in the R k are less than the 
uncertainties in the prior. The ratio Ay/Aj = 1.63 ± 0.05 lies between the Key project 
{AyjAi = 1.6) and SN la {AyjA 1 = 1.7) values. The shifts in the Glass-Carter JHK 
coefficients were expected, because we lacked the true mean effective wavelengths for these 
filters and had simply set their values to those for the CTIO JHK filters. However, the 
sign of the shift is opposite to that found by Laney & Stobie (1993), who found smaller 
extinction coefficients using a different analysis technique. 

Our distance estimates generally agree with the published values (see Tables 3 and 4, 
and Figures 1 and 2a) with a few exceptions. We find that M 31, M 33, and NGC 300 are 
closer than the published values by —0.11, —0.12, and —0.05 magnitudes respectively. In 
M 33 the reason is the significantly higher extinction estimate of E(B — V) — 0.16 instead 
of 0.10. Both earlier Cepheid distance estimates for M 33 (Freedman 1985, Madore et al. 
1985), and other models tested in Freedman et al. (1991) match our higher extinction 

1 Laney & Stobie (1994) found zero-points of V ok = 13.28 + 0.09, 11.90 + 0.06, 11.47+0.06, 
and 11.38 ± 0.06 mag and PL slopes of a' k = -2.87 ± 0.07, -3.31 ± 0.05, -3.42 ± 0.05, and 
—3.44 ± 0.05 in the V, J, H, and K bands for a mean LMC modulus of 18.5 mag 

2 Tanvir (1996) found zero-points of V ok = 13.24 ± 0.04 and 12.45 ± 0.03 and PL slopes 
of a' k = -2.774 ± 0.083 and -3.039 ± 0.059 for V and I. 
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estimates. Since the Cepheid distances to M 31 and M 33 calibrate many other distance 
indicators, these shifts would increase some estimates of the Hubble constant by 6-7%. 
Our agreement with the Key Project estimates is excellent, with a mean shift of —0.05 
mag, similar to that found by Tanvir (1996) due to zero-point recalibration. Our error 
and extinction estimates are also in good agreement, except for M 81 where we find a 
significantly higher mean extinction of 0.075 (the Key Project value of 0.03 is less than 
the estimated foreground extinction of 0.04). For M 96, we agree with the values for the 
distance and mean extinction found by Tanvir et al. (1995), but we derive much larger 
uncertainties of 0.28 (versus 0.16) mag for the distance and 0.12 (versus 0.03) for the 
extinction. We found significant disagreements with the Type la project in the distance 
estimates, extinction estimates and uncertainties. On average our distances were —0.09 
mag smaller, the extinctions were 0.04 larger, and our distance uncertainties were generally 
twice as large. The results significantly change the zero-point of the MLCS SNIa distance 
scale (Riess et al. 1996), since the distances to the calibrators 1972E in NGC 5253, 198 IB 
in NGC 4536, and 1990N in NGC 4639 are revised from 28.08 ± 0.10 mag to 27.70 ± 0.32, 
31.10 ± 0.13 to 30.97 ± 0.22, and 32.03 ± 0.22 to 32.11 ± 0.32 respectively. The revised 
MLCS estimate of the Hubble constant is 69 ± 8 km s" 1 Mpc -1 , an 8% increase from the 
original estimate of 64 ± 6 km s" 1 Mpc^ 1 . Some of the differences may be due to our use 
of a consistent extinction law for all the galaxies (the Type la project uses the Madore & 
Freedman (1991) PL relation based on Ry/Ri = 1.6 but fit the data using Ry/Rj = 1.7), 
and our use of statistically consistent treatments of the covariance between distance and 
extinction. For example, in NGC 5253 Saha et al. (1995) give large uncertainties in the 
extinction (see Table 3) that are statistically inconsistent with the small uncertainties in the 
distance. Figure 2 compares two absolute distance estimates, the expanding photosphere 
method for Type II supernovae (Eastman et al. 1996) and physical models of Type 
la supernovae (Hoflich & Khokhlov 1996), and two relative distance estimates, surface 
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brightness fluctuations (Tonry 1996) and the MLCS/SNIa (Riess et al. 1996) method, to 
the Cepheid distances (see Table 4). 

Magnitude selection biases may be important in Model because of the large scatter 
in the PL relations. Biases have been extensively discussed in the literature, most recently 
by Tanvir (1996), along with debates on how to optimally perform model fits (e.g. fitting 
inverse PL relations as in Kelson et al. (1996)). We can avoid the biases almost completely 
by using PLC relations to model the correlations in the residuals and to reduce the intrinsic 
scatter. Since these approaches are also physically more useful, we will not discuss selection 
biases further. 

3.2. Model 1: Scatter In the Extinction 

The magnitude residuals of all the extragalactic Cepheids are highly correlated and 
lie along the extinction direction (see Figure 3).0 The importance of allowing individual 
extinctions is glaringly obvious if we compare the residuals for the Galactic and extragalactic 
samples - for the Galactic Cepheids there are 624 measured mean magnitudes for 121 
Cepheids with an rms residual of 0.06 magnitudes, while for the extragalactic Cepheids 

3 Contrary to the conclusions of the Type la project, we find that the residuals in IC 4182, 
NGC 4536, NGC 4496a and NGC 5253 are also correlated with the extinction vector. We 
believe their conclusion is an artifact of examining residuals in the space of AV — Al versus 
AV. Both the mathematics of the covariance matrix and simple experiments show that 
when the uncertainties in / are larger than those in V, the AV — Al versus AV covariance 
diagrams can look symmetric and uncorrelated even when the data is correlated. The 
quantity AV —Al simultaneously reduces the signal from differential extinction and increases 
the apparent noise (see Fig. 3). 
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Fig. 1. — Distance comparisons between Model and either published (top) or Model 1 
(bottom) distances. The horizontal error bar is the error in Model 0, while the vertical error 
bar is the error in the comparison distance. 
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Fig. 2. — Distance comparisons to other methods. The expanding photosphere method 
for SNII (EPM, solid squares, Eastman et al. 1996) and physical models of SNIa (Snla, 
solid triangles, Hoflich & Khokhlov 1996) are absolute distances. The surface brightness 
fluctuation distances are shown relative to M 31 (SBF, open squares, Tonry et al. 1996) and 
the SNIa/MLCS method distances are shown relative to NGC 4536 (MLCS, open triangles, 
Riess et al. 1996). The error bars are the uncertainties in the non-Cepheid distance indicator. 
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Fig. 3. — Magnitude residuals in Model 0. The left panel is a simple diagram of the AV and 
AI residuals (squares). It shows that the residuals are strongly correlated and lie roughly 
along the extinction vector. The right panel displays the same residuals using the AV 
versus AV-AI format of Type la project. The correlation, while still clearly present, is 
much more difficult to recognize. We also display the residuals (vertically offset, crosses) 
and the extinction direction for the more accurately measured LMC and SMC Cepheids to 
illustrate the effects of measurement uncertainties on the apparent correlations. The heavy 
solid, dashed and light solid lines show the expected directions for distance, temperature and 
extinction residuals respectively. 
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there are 1617 magnitudes for 545 Cepheids with an rms residual of 0.29 magnitudes. 
The independent distance and extinction estimates for the Galactic Cepheids dramatically 
reduce their residuals compared to the rest of the sample even though a large subset of the 
Galactic Cepheids have 5 to 8 color photometry while most of the extragalactic Cepheids 
have only two color photometry. 

In Model 1 we allow all Cepheids an independent extinction variable, as we have 
already used for the Galactic Cepheids in Model 0. Each Cepheid % in galaxy k was assigned 
extinction (E^) + AE^, where (E^) is the mean extinction and the correction AE^ is forced 
to have zero mean for each galaxy and is limited by a Gaussian prior whose width OEk is 
simultaneously optimized. We need the prior because most of the extragalactic Cepheids 
have only two filters, and as we start assigning each Cepheid individual extinctions, 
temperatures and metallicities we will overfit the data if we do not include the covariance 
matrix (eqns. 4 and 5) and optimize the priors. Some extragalactic distance estimates have 
used of individual extinction estimates (e.g. Tanvir et al. 1995) or reddening free magnitude 
estimates (e.g. Freedman et al. 1990, 1991, 1992, Saha et al. 1996ab, 1997), although the 
final distance estimates are usually based on the standard PL relations. These individual 
extinction estimates correspond to taking the limit oe — ► oo in eqns. (4) and (5), and they 
may overestimate the width of the extinction distribution in noisy data. 

The changes in the Cepheid relations and distances are summarized in Tables 2 and 3 
and Figure 1. In Model 1 the rms residuals for the extragalactic Cepheids drop to 0.09 mag 
from 0.29 mag in Model 0. The rms residuals for the LMC and SMC (0.06 and 0.08 mag) 
are smaller than the other extragalactic systems (0.10 mag) even though the average Cloud 
Cepheid was measured in 5 filters. The removal of the extinction-correlated residuals leads 
to large reductions in the PLC relation widths (which represent uncorrelated errors), but 
the residuals are still correlated (particularly UJHK) and the PLC widths are still broader 



-27- 



than expected given the estimated measurement errors. The zero-points and slopes show 
little change from Model and have reduced uncertainties, but the estimated extinction 
vector has changed considerably in the infrared, suggesting that it has absorbed some of 
the variance due to temperature as well as extinction. In particular, the distance of the 
sun from the Galactic center (R ) is strongly affected by the change in the IR extinction 
coefficients (from 7.5 kpc to 6.7 kpc). 

As Figure 3 illustrates, the extinction and temperature are very nearly parallel for V 
and I, so it is very difficult to cleanly separate the two physical terms in noisy two-color data. 
Madore & Freedman (1991) question whether the two variables can be accurately separated 
even with more accurate three-color data used for Galactic Cepheid extinction estimates, 
although their position is strongly rejected by Laney & Stobie (1993, 1994). Moreover, the 
HST Cepheid data may have correlated systematic errors due to crowding, differences in 
analysis methods, and the construction of the I mean magnitude using the V light curves 
for interpolation that Model 1 interprets as extinction variations. Such systematic errors 
have no effect on the distance estimates, since correlated residuals must be modeled for a 
correct statistical treatment, but they may strongly affect our interpretation of the scatter 
in terms of extinction. Note, however, that the spread in color in the extragalactic systems 
is quite comparable to the spread in the Magellanic Clouds (Figure 3), and it would be 
truly astonishing if the extinction from a large sample of objects randomly chosen from 
spiral galaxies failed to show significant scatter. 

Other than the value of Ro, the distances are little changed from Model (see Fig. 
1), with a mean shift of less than 0.01 mag. The typical distance and mean extinction 
uncertainties, however, are roughly half those of Model 0, although the reduction is most 
dramatic for galaxies with small numbers of Cepheids (e.g. M 96). The change is largely 
due to the reduced rms magnitude residual in estimating the intrinsic Cepheid luminosities. 
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In Model the statistical uncertainty is ~ 0.29/iV 1 / 2 magnitudes where the numerical 
coefficient is the combination of the measurement errors and the PLC relation width, while 
in Model 1 the typical uncertainty is ~ 0.10/iV 1 / 2 magnitudes. The total uncertainties 
are larger due to calibration uncertainties in the PLC relations and the HST magnitudes 
that are unaffected by statistical averaging. Fitting the correlated residuals also makes our 
model significantly less sensitive to magnitude selection biases than Model 0, although the 
general agreement of the results suggest they were not of great importance. The MLCS 
estimate of the Hubble constant becomes 72 ± 6 km s _1 Mpc -1 , slightly higher than in 
Model 0. Figure 2 compares the Cepheid distances to the distances in Table 4. 



3.3. Model 2: Temperature and the Positivity of the Extinction 

Figure 4 shows the mean extinctions relative to the foreground extinction and the 
width of the extinction distribution as a function of the metallicity of the host galaxy. The 
Model 1 solutions are unphysical because many of the galaxies have either negative mean 
internal extinctions (SMC, NGC 300, IC 4182 and NGC 4639) or extinction distributions 
extending to negative internal extinctions. Even for the LMC and the SMC, where we 
have many filters, good accuracy, and no systematic problems such as the HST calibration 
terms, there are Cepheids with extinctions less than the expected foreground extinction. 
The simplest solution to the problem, and the one advocated by Freedman et al. (1992) 
when they first faced the problem for the NGC 300 Cepheids and more recently advocated 
by Bohm-Vitense (1997), is simply to raise the mean LMC extinction from the standard 
{E) LMC = 0.10 assumed in the extragalactic Cepheid analyses. Both Freedman et al. 
(1992) and Bohm-Vitense (1997) suggest that the problem can be solved by increasing 
the mean extinction to (E)lmc — 0.18. Our larger sample and the interpretation of the 
scatter in color as an extinction distribution increases the required LMC extinction to 
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(E) LMC 0.25. The standard Galactic Cepheid analyses derive mean LMC extinctions 
of only 0.07 (e.g. Caldwell & Coulson 1986), even lower than {E) LMC = 0.10, based on a 
combination of the Cepheid colors and other estimates of the space reddening. Grieve & 
Madore (1986) found median LMC supergiant extinctions of approximately 0.10, with 90% 
of the estimated extinctions below E(B — V) = 0.18. Bessel (1991), in arguing for higher 
than generally accepted extinctions for the LMC and SMC, advocates an increase to only 
(E) ~ 0.13 based on HI column densities, optical polarization, interstellar absorption lines, 
and color excesses. Thus the trivial solution to the negative extinction problem of simply 
raising {E) LM c can only be a partial solution. As Freedman et al. (1992) also noted, the 
temperature distribution at fixed period contributes to the negative extinction problem if 
the two variables are not distinguishable and hotter than average Cepheids are assigned 
reduced extinctions. 

In Model 2 we try to minimize the effects of the temperature distribution in biasing 
the extinctions by giving each Cepheid a temperature deviating by St from the mean for its 
period, constrained by a Gaussian prior of width <jpc corresponding to the width of the 
instability strip (see eqn. 2). We also force the extinctions to be positive by adding the 
extra terms to the likelihood function discussed in §2.2. Since there is some uncertainty 
in the extinction normalization, we explore three models with {E) LMC = 0.10, 0.15, and 
0.20 labeled by Model 2-10, 2-15, and 2-20, that correspond to low, slightly high, and very 
high estimates for the mean extinction. We adopt Model 2-15 as our standard, and the 
resulting PLC vectors, distances, and extinctions are presented in Tables 2 and 3. In Model 
2 we must determine the temperature vector (3k as well as the temperature corrections, St, 
and as is frequently noted in the extragalactic critiques of the Galactic Cepheid methods 
(e.g. Madore & Freedman 1991), the temperature state vector (3k is hard to determine 
uniquely from the mean magnitudes. Indeed, we found that we could not stably estimate 
(3k without including a prior based on the variation of color with phase. Based on an 
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analysis of Cepheid light curves (Kochanek 1997b), we estimated a model (3 0k (see Table 2), 
assumed that the uncertainty in the coefficients was 0.05, and added a Gaussian prior for 
the deviation of f3 k from p ok to the likelihood. We should note, however, that many of the 
Galactic analyses use a very similar approach to calibrating the PLC relations (e.g. Laney 
& Stobie 1986, 1993, 1994). The need to make such a strong assumption about f3 k combined 
with (or due to) the limited two-color photometry on most of the extragalactic systems is a 
primary limitation for the remainder of our analysis, but the simultaneous optimization of 
the prior widths keeps us from overfitting the data. 

We find a width for the instability strip of a PC ^ 0.18 ± 0.01 in all three models, 
which corresponds to a FWHM in the V-I (B-I) colors at fixed period of 0.08 (0.14) mag. 
Unlike the previous two models, Model 2 gives non-zero values for HST calibration variables 
(AV,AI) of (0.05±0.04,-0.04±0.04), (0.03±0.04,-0.03±0.04), (0.02±0.04,-0.03±0.04) 
for Models 2-10, 2-15, and 2-20 respectively. In all three cases the calibration uncertainties 
provide a significant part of the solution to the negative extinction problem. Note that 
the PLC zero points for Model 2-15 are not directly comparable to the previous models 
due to the change in the extinction. The infrared extinction coefficients are now closer to 
the priors than in Model 1 and agree with the Laney & Stobie (1993) estimates and the 
estimated value of R = 7.7 ± 0.3 kpc is again consistent with other estimates (Reid 1993). 

Standard Cepheid distances are biased distance estimators because they treat positive 
and negative extinctions equally. The sense of the bias is always to overestimate the 
distances, because lower extinctions correspond to higher distances. Random photometry 
errors can and will produce negative extinctions for objects with positive extinctions, 
but such problems affect only the mathematics of implementing the positive extinction 
condition, not the need for it. Any procedure to impose the physical condition that the 
extinction is positive must reduce the likelihood for large distances and thereby drive the 
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Fig. 4. — Extinctions and metallicity. The points for each galaxy show the mean extinction 
relative to the foreground extinction (Ek) — Efk and the width of the extinction distribution 
OEk for Model 1 as a function of the metallicity. The uncertainty in the mean extinction is 
smaller than a eh by roughly the square root of the number of Cepheids in the galaxy (see 
Table 3). The metallicity error bars show the spread in the metallicities of the Cepheids used 
in the model. The solid and dashed lines show the best linear and constant fits for the trend 
excluding NGC 300 and NGC 4639. Several points have been shifted slightly in metallicity 
to separate the vertical error bars. 
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Fig. 6. — Metallicity corrections to the zero point. Likelihood contours for the change in 
mean magnitude ((7y + 7j)/2) and V-I color (7(7 — 7/) compared to other observational and 
theoretical determinations. The contours are the 1-a (A% 2 = 2.30) and 2-<r (A% 2 = 6.17) 
confidence intervals for two parameters. The G94 (Gould 1994) and MF90 (Madore & 
Freedman 1990, as estimated by Gould 1994), estimates included no color variation. CWC93 
marks the theoretical estimate for V-I from Chiosi et al. (1993). S88 and St95 mark the 
theoretical estimates for B-V from Stothers (1988) and Stift (1990, 1995). CC87 marks the 
semi-empirical model for B-V from Caldwell & Coulson (1987). S96 marks the Sasselov et 
al. (1996) V and I determination from the EROS Cepheids. 
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best fit distance estimate downwards. It is also a global bias, because the Cepheid distance 
and extinction estimates are tightly correlated by the PLC relations. When we force a 
galaxy with strongly negative extinction estimates like NGC 300 to have a higher extinction 
and a lower distance, the correlations caused by all the Cepheids sharing the same estimate 
of their intrinsic colors will also force galaxies like M 33 with strictly positive extinctions 
to higher extinctions and lower distances as well. Adjusting the extinction of one galaxy 
without adjusting the extinction of all galaxies is physically and statistically incorrect, unless 
the reason for the adjustment is a systematic error in the data for that particular galaxy. 
The negative extinction problems should not be solved "locally." 

As Figure 5 shows, Model 2 simply drives the galaxies to lower distances and higher 
extinctions, until it is reasonably unlikely that any Cepheids have negative extinctions. 
The mean distances shift by —0.20, —0.14, and —0.05 mag for Models 2-10, 2-15, and 
2-20 respectively, with the low LMC extinction Model 2-10 requiring the largest shift. The 
effect is strongest on the low metallicity galaxies, which means the changes in the MLCS 
calibrations are considerable. The MLCS Hubble constant estimates are 80 ±6, 78 ± 6, 
and 72 ± 6 km s _1 Mpc -1 for the three models compared to 72 ± 6 km s _1 Mpc -1 for 
Model 1. Because the effects of the bias are global, even a relatively high extinction galaxy 
like M 100 is driven to lower distances. Scaling from the Mould et al. (1995) estimate of 
80 km s -1 Mpc -1 , which becomes 82 km s -1 Mpc -1 in Model 1, the Hubble constant shifts 
to H = 90, 88, and 84 km s" 1 Mpc" 1 in Models 2-10, 2-15, and 2-20 respectively^ Despite 
the large shifts, the agreement with other distance indicators (see Fig. 2) is no worse than 
in Model 1, if only because the absolute distance indicators have low accuracy while the 
relative distances are almost unchanged. The goodnesses of fit for Model 2 are considerably 

4 The uncertainties in the H estimate from M 100 are dominated by the systematic 
uncertainties in the model for Virgo (Mould et al. 1995). 
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worse than Model 1, even though we have allowed separate temperature estimates for the 
Cepheids. The likelihood function has increased by A21nL = 139, 78, and 65 relative to 
Model 1, despite adding 695 new degrees of freedom to the model! 

3.4. Model 3: Metallicity 

The Model 2 solutions are not very attractive because of the huge downward shifts in 
the distances, and the association of unphysical extinctions with low metallicity in Figure 4 
strongly suggests the need for a metallicity dependence in the Cepheid distance. In Model 
3 we add a metallicity dependence to the Cepheid zero-point vector j' k , as in the model of 
Sasselov et al. (1996). The formal trend of extinction with metallicity (neglecting NGC 300 
and NGC 4639) corresponds to a color dependence of V-I oc (0.09±0.01)A[O/#], which is 
comparable to the color variations predicted in theoretical models (e.g. Stothers 1988, Stift 
1990, Chiosi et al. 1983) and previous experimental estimates (e.g. Gieren et al. (1993), 
Sasselov et al. (1996)). There are, however, two discrepant points in Figure 4, NGC 300 
and NGC 4639, both of which are significantly bluer than the trend. For NGC 4639 the 
large uncertainties in the extinction can probably explain the discrepancies. NGC 300 is 
peculiar, and we drop it from our subsequent analysis - in doing so, we are assuming that 
the blue color of the Cepheids is either an artifact of the Freedman et al. (1992) data or 
that the HII region abundance estimates are incorrect. 

In Model 3 we added the metallicity dependent zero-point correction 7^AZ to Model 
2, and repeated the calculation for the three different mean LMC extinctions labeled by 
Model 3-10, 3-15, and 3-20. As in Model 2, we present the full results only for Model 3-15 
in Tables 2 and 3. We fit the data leaving the input metallicities fixed (a z — > 0), because 
attempting to determine the metallicities from the data was unjustified by the amount and 
quality of the data. If we have found a solution for the distance modulus ji\ and extinction 
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Ei ignoring the metallicity zero point correction 7^, then the true distance modulus is 



= |Xl -T^AZ+^: y ^ y + ^AZ (9) 

and the true extinction is 



E = E 1 - l v ^ AZ (10) 

where the terms are broken into the change in the mean V and I magnitudes and the 
change in the V-I color. If low metallicity Cepheids are blue — 7/ > 0) then adding 
the metallicity term can solve the extinction problem without driving the distances of all 
galaxies downwards as in Model 2. The metallicity terms will not, however, change the 
embarrassingly low distances to the low metallicity galaxies in Model 2, because their 
distances are still reduced by (R v + -R/)/2 ~ 2.7 times the change in the extinction. 

Our solutions generically made the metal rich Cepheids redder and brighter than the 
metal poor Cepheids, as shown in Figure 6. The best fit solution changes little with the 
assumptions about the LMC extinction, and the absence of a metallicity dependence is 
ruled out at greater than 95% confidence. The mean luminosity changes in the V and I 
bands, {i v + 7^), are -0.13 ± 0.09, -0.15 ± 0.14, and -0.17 ± 0.15 mag/dex for Models 
3-10, 3-15, and 3-20 respectively, and the mean changes in color, 7^ — 7^, are 0.15 ± 0.03, 
0.13 ± 0.04, and 0.15 ± 0.07 respectively. The less strongly we restricted the permitted 
range for the extinctions, the more uncertain the metallicity dependence. As pointed out 
by Sasselov et al. (1996), the positivity of the extinction is an important component in 
the quantitative determination of composition effects. The HST calibration uncertainties 
still represent a major problem, and they contribute much of the solution to maintaining 
positive extinctions. The calibration variables (AV, AI) are (0.09 ± 0.04, -0.07 ±0.04), 
(0.05 ± 0.04, -0.05 ± 0.04), and (0.03 ± 0.04, -0.03 ± 0.04) for the three solutions. The 
overall structure of the metallicity vector 7^ matches theoretical expectations. Metal rich 
Cepheids show a decreased flux in U and B, and then a gradually increasing flux in the 
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redder bands, reaching a plateau in the infrared. The uncertainties in the 7^ in Table 2 are 
dominated by the uncertainties of the mean luminosity change, and as the error ellipses in 
Figure 6 demonstrate, the uncertainties in the color changes are considerably smaller. 

Most of the debate about composition effects on extragalactic Cepheid distances has 
focused on the M31 Cepheids studied by Madore & Freedman (1990) and Gould (1994). 
These models allowed no color changes due to composition (7^ = 7^), while it is clear 
both from theoretical models, previous experimental estimates, and our estimates, that the 
changes in color are the dominant source of changes in distance estimates. Our estimates of 
the effect are slightly higher than the theoretical estimates of Chiosi et al. (1993). The B-V 
color change of 0.28 mag/dex is not as well constrained because most of the galaxies lack 
B photometry, but it is larger than the theoretical estimates of Stothers (1988) and Stift 
(1990, 1995). The values match the experimental determination by Sasselov et al. (1996) 
using the EROS sample of Cepheids in the LMC and SMC and the typical Galactic Cepheid 
metallicity correction models (e.g. Caldwell & Coulson 1986, 1987, Gieren et al. 1993). 

Figure 7 shows the changes in the distances relative to Model 1 for the three different 
assumptions about the LMC extinction. As expected, the metal poor galaxies (e.g. 
NGC 5253, IC 4182, NGC 4536) are shifted to lower distances relative to the metal rich 
galaxies (e.g. NGC 3351, M 33, M 100). The mean change in distance depends on the 
mean extinction of the LMC, with Model 3-10 still requiring a significant reduction in the 
mean distances. The MLCS estimates of the Hubble constant become 85 ± 6, 80 ± 6 and 
79 ± 6 km s" 1 Mpc" 1 for Models 3-10, 3-15, and 3-20 respectively, while the estimates 
based on M 100 (Mould et al. 1995) are 83, 78 and 74 km s" 1 Mpc" 1 respectively. The 
metal poor calibration (MLCS) moves to lower distance and higher Hubble constant, while 
the metal rich calibration (M 100) moves to higher distances and lower Hubble constants. 
The comparison with other distance indicators (Figure 2) is no worse than the other cases, 
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with the exception of NGC 5253. 

4. Conclusions 

We have systematically explored the problem of extragalactic Cepheid distances and 
their primary systematic errors. While the details of some of the model implementations 
are certainly open to criticism, we have tried to include or illustrate all the principal 
uncertainties. When we implement the standard extragalactic analysis method (Madore 
& Freedman 1991) in Model 0, we find general agreement with existing distances, with 
significant corrections only for the Type la supernova calibration galaxies. Our revised, 
distances to the SNIa MLCS (Riess et al. 1996) calibrating galaxies NGC 5253, NGC 4536, 
and NGC 4639 produce a revised Hubble constant estimate of 69 ± 8 km s" 1 Mpc" 1 . Our 
new distance estimates and uncertainties have the advantages of a homogeneous statistical 
treatment, a standard extinction model, and the inclusion of the full uncertainties in the 
extinction and Cepheid model on the distances. In Model 0, the magnitude residuals for all 
galaxies are strongly correlated with the extinction vector, and when we allow each Cepheid 
an individual extinction in Model 1, the typical magnitude residual drops from 0.29 mag to 
0.09 mag. The Galactic Cepheid community (e.g. Caldwell & Coulson 1986, 1987, Caldwell 
& Laney 1991, Fernie 1990, Fernie et al. 1995, Laney & Stobie 1993, 1994, 1996) includes 
individual extinctions as a matter of routine. Some extragalactic studies effectively fit 
individual extinctions or use reddening free magnitudes (e.g. Tanvir et al. 1995, Freedman 
et al. 1990, 1991, 1992, Saha et al. 1996ab, 1997), but they are usually not used in the 
final distance estimate. As the sadly neglected work of Gould (1994) emphasized, a correct 
statistical model must include the effects of these correlations on the uncertainties in the 
distance independent of the physical interpretation for their origin. The danger of confusing 
extinction, temperature, and correlated systematic errors affects only the interpretation 
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of the model not the need to account for correlated residuals. Including the correlations 
does not change the values for the distances and mean extinctions, but it does significantly 
reduce their uncertainties (see Figure 1). 

In both of these models, many Cepheids require negative intrinsic extinctions for the 
host galaxies, which is clearly unphysical. The simple solution to the problem is to raise 
the mean LMC extinction form the standard value of (E) lmc = 0.1 to ^ 0.2, as first 
suggested by Freedman et al. (1992) to solve the negative extinction estimates in NGC 300 
and recently reintroduced by Bohm-Vitense (1997) in an analysis of the Galactic, LMC, 
SMC, and M31 Cepheids. In our larger sample and after including the scatter in the 
extinction, we find the problem is significantly worse and that a mean LMC extinction of 
^ 0.25 would be required to eliminate the problem completely. However, even the high 
estimates of the mean extinction in the LMC by Bessel (1991) based on polarization, HI 
column density, interstellar absorption lines, and intrinsic color estimates correspond only 
to (E) L mc ^ 0.13, and 90% of the LMC supergiants studied by Grieve & Madore (1986) 
had extinctions less than 0.18, so simply raising the mean LMC extinction is an implausible 
solution to the problem of negative extinctions. If we can reject raising the LMC extinction, 
then the need for negative extinctions is conclusive evidence for adding additional physics 
to the Cepheid model or for substantial correlated systematic errors in the magnitudes 
beyond the simple HST calibration uncertainties. 

The physical requirement that the extinction be positive also means that the standard 
Cepheid distance estimates are systematically biased by their equal treatment of positive 
and negative extinctions. The covariance of distance and extinction mean that low 
extinction estimates are associated with higher distances, so any requirement for positivity 
in the extinction will drive the distances to all Cepheid galaxies downwards. The bias 
extends equally to galaxies with and without a negative extinction problem because the 
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extinction estimates and distances are also correlated between galaxies - if I force a 
reduction in the distance to one galaxy, the correlations lead to a reduction in the distances 
to all galaxies. 

A partial solution to the negative extinction problem is to use a finite temperature 
distribution at fixed period so that blue Cepheids are interpreted as being hotter rather 
than having lower extinctions (see the discussion in Freedman et al. 1992). It is only 
a partial solution both because the instability strip is narrow and because it cannot 
significantly change the mean extinction of the distribution unless there are few Cepheids 
in the galaxy. In Model 2 we allowed the Cepheids a temperature distribution, and forced 
the positivity of the extinction for a range of LMC mean extinctions. As expected from the 
bias in the standard distance estimates, the distance to every Cepheid galaxy decreased by 
mean values of —0.20, —0.14, and —0.05 mag for LMC mean extinctions of {E) LMC = 0.10, 
0.15, and 0.20 respectively. All Hubble constant estimates systematically rise, although the 
effect is more dramatic for the low metallicity Type la Project galaxies because of their 
bluer Cepheids. In particular, the MLCS Hubble constant estimates become 80 ± 6, 78 ± 6 
and 72 ± 6 km s' 1 Mpc" 1 and the M100 (Mould et al. 1995) estimates become 90, 88, and 
84 km s^ 1 Mpc -1 in order of increasing LMC extinction. Model 2 is only a partial solution 
to the negative extinction problem, and despite the addition of a temperature variable for 
all 694 Cepheids, the likelihood of the solutions declined significantly compared to Model 1. 

The extinctions shows a rough correlation with the metallicity of the host galaxy, as 
would be expected from theoretical models predicting that metal poor Cepheids should 
be bluer than metal rich Cepheids (e.g. Stothers (1988), Stift (1990, 1995), Chiosi et al. 
1993). In Model 3 we added a metallicity correction to the Cepheid zero-point, similar to 
the model of Sasselov et al. (1996), and we find a change in the mean V and I magnitude of 
— 0.14±0.14 mag/dex and a change in the V-I color of 0.13±0.04 mag/dex. The metallicity 
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terms change little with the assumptions about the LMC extinction. The trend with 
wavelength is that metal rich Cepheids become fainter in U and B, start getting brighter 
at V, and show the largest increases in the infrared, as expected from line-blanketing 
increasing the opacity towards the blue and back-warming increasing the emission towards 
the red. While the quantitative estimates may be somewhat high, they roughly agree with 
both the theoretical expectations and other experimental or semi-empirical determinations 
(e.g. Caldwell & Coulson 1986, 1987, Gieren et al. 1993, Stift 1995, Sasselov et al. 1996). 
The change in color with metallicity is more important than the change in luminosity, so 
studies focusing only on the change in luminosity (Madore & Freedman 1990, Gould 1994) 
miss the dominant effect. The addition of the metallicity correction lets the model solve 
the negative extinction problems without simply driving the distance estimates downwards, 
and as discussed in Gould (1994) and Sasselov et al. (1996), it can also explain many of the 
discrepancies between low (Type la supernovae in low metallicity galaxies) and high (other 
distance indicators in high metallicity galaxies) estimates of the Hubble constant. The 
MLCS Hubble constant estimates from the low metallicity galaxies are 85 ± 6, 80 ± 6 and 
76 ± 6 km s" 1 Mpc^ 1 , while the high metallicity M100 estimates are 83, 78 and 74 km s^ 1 
Mpc" 1 as we increase the mean LMC extinction. A metallicity dependence to the Cepheid 
distance scale is the most plausible explanation of the negative extinction problem if it 
cannot be explained by systematic errors in the photometry. Otherwise, the qualitative 
effect is probably secure even if our quantitative results are changed by improved data or 
metallicity estimates 

There is no question that the models including temperature and metallicity are 
beginning to push the limits of the data, although our inclusion of the covariances 
between variables should properly treat near degeneracies. We can identify five areas for 
improvement in the data. The first is that the data on the LMC Cepheids needs to be 
expanded and observed in a more uniform set of filters. The enormous MACHO and EROS 
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samples of Cepheids can almost fill this role, but the calibration uncertainties need to be 
reduced, comparisons to standard filters need to be better understood, and comparable 
data in additional filters is required. The second problem is that the HST Cepheid samples 
have only two-color photometry. Since even the simplest models that can fully probe the 
systematic uncertainties in the Cepheid distance scale require extinction, temperature and 
composition estimates for each individual Cepheid, a minimum of four-color photometry 
is needed. The two additional filters should be B (sensitive to extinction and metallicity) 
and H or K (sensitive to temperature and metallicity). Without additional colors it will be 
difficult or impossible to cleanly address the systematic problems in the Cepheid distance 
scale. The third problem is that the formal uncertainties in the absolute calibration of 
the HST magnitudes (about 0.05 mag) are large compared to the effects we are probing. 
Some means of improving the absolute calibrations is critical to better controlling the 
systematic uncertainties in the Cepheid distance scale. The fourth problem is that many 
of the Cepheid host galaxies lack metallicity measurements, forcing us to rely on general 
correlations between galaxy type or luminosity and metallicity. 

The final problem is that no matter how elaborate or rococo we make the analysis 
procedures for Cepheid mean magnitudes, it is not the proper way to analyze the data. 
The only "correct" way to analyze the Cepheid data is to directly fit the observed light 
curves, in the spirit of the MLCS method for Type la supernovae (Riess et al. 1995). Such 
an approach is critical to the extragalactic samples where heroic investments of HST time 
still produce poor phase coverage in V and terrible phase coverage in I, let alone adding 
data in the two additional filters needed to control systematic uncertainties. Stetson (1996) 
has developed a template method for identifying Cepheids and better estimating mean 
magnitudes, and in our second paper on Cepheid distances (Kochanek 1997) we develop 
an improved template method incorporating both color and velocity data. Ideally, the 
shapes of the light curves also constrain the gross physical properties of the Cepheids and 
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break some of the close degeneracies seen between temperature, extinction, distance, and 
metallicity that limit the validity and accuracy of our current conclusions. 
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Table 1. Cepheid Data 
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Note. - The metallicities are from Zaritsky et al. (1994), except for M 96 (Oey & Kennicutt 
1993) and M 31 (Blair et al. 1982). Metallicities for the galaxies not included in Zaritsky et al. 
(1994) were estimated from the metallicity-type or metallicity-magnitude relations. The metallicity 
uncertainties represent the rms range of metallicities assigned to the Cepheids in that galaxy based 
on the metallicity gradients if known. The foreground extinction estimates Ef = E(B — V)f are 
from Burstein k, Heiles (1984). 



Table 2. Cepheid Correlations 
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f 

a k 




— 2.43 ± 0.14 


— 2.76 ± 0.11 


— 2.94 ± 0.09 


— 3.06 ± 0.07 












VPLCk 




0.36 


0.27 


0.22 


0.18 








MF2 




Vk 




14.03 ± 0.16 


13.23 ± 0.12 


12.80 ± 0.11 


12.40 ± 0.09 


11.89 ± 0.07 


11.50 ± 0.06 


11.43 ± 0.05 






1 

a k 




— 2.53 ± 0.28 


— 2.88 ± 0.20 


— 3.04 ± 0.17 


— 3.14 ± 0.17 


— 3.31 ± 0.11 


— 3.37 ± 0.10 


— 3.42 ± 0.09 






CPLCk 




0.40 


0.29 


0.25 


0.21 


0.16 


0.14 


0.13 





4277 


v k 


14.54 ± 0.05 


14.05 ± 0.03 


13.28 ± 0.03 


12.83 ± 0.03 


12.45 ± 0.03 


11.89 ± 0.02 


11.53 ± 0.02 


11.45 ± 0.02 






1 

a k 


O C\i~' 1 C\ "1 A 

— 2.06 ± 0.14 


— 2.30 ± 0.08 


— 2.62 ± 0.05 


— 2.81 ± 0.06 


— 2.95 ± 0.05 


— 3.22 ± 0.05 


— 3.31 ± 0.05 


— 3.35 ± 0.04 






Rk 


5.20 ± 0.08 


4.30 ± 0.06 


= 3.3 


2.66 ± 0.05 


2.02 ± 0.05 


1.05 ± 0.06 


0.66 ± 0.06 


0.47 ± 0.06 






CPLCk 


0.27 ± 0.02 


0.29 ± 0.01 


0.27 ± 0.01 


0.15 ± 0.01 


0.21 ± 0.01 


0.13 ± 0.01 


0.11 ± 0.01 


0.10 ± 0.01 


1 


7418 


v k 


14.60 ± 0.02 


14.06 ± 0.01 


13.29 ± 0.01 


12.86 ± 0.01 


12.44 ± 0.01 


11.86 ± 0.01 


11.49 ± 0.01 


11.41 ± 0.01 






f 

a k 


— 1.83 ± 0.12 


— 2.25 ± 0.05 


— 2.63 ± 0.07 


— 2.79 ± 0.03 


— 2.95 ± 0.04 


— 3.16 ± 0.03 


— 3.26 ± 0.03 


— 3.30 ± 0.03 






Rk 


5.13 ± 0.04 


4.28 ± 0.02 


= 3.3 


2.72 ± 0.01 


2.17 ± 0.02 


1.31 ± 0.03 


0.95 ± 0.03 


0.78 ± 0.03 






CPLCk 


0.03 ± 0.03 


0.00 ± 0.01 


/~\ r\r\ 1 /~\ m 

0.00 ± 0.01 


0.00 ± 0.01 


0.02 ± 0.01 


/ _ \ /^/^ 1 r\ r\-i 

0.06 ± 0.01 


0.06 ± 0.01 


0.07 ± 0.01 


2-15 


7340 


v k 


14.34 ± 0.03 


13.82 ± 0.02 


13.13 ± 0.02 


12.75 ± 0.02 


12.38 ± 0.02 


11.87 ± 0.02 


11.53 ± 0.02 


11.46 ± 0.02 






1 

a k 


— 1.95 ± 0.12 


— 2.31 ± 0.09 


— 2.67 ± 0.07 


— 2.80 ± 0.06 


— 2.96 ± 0.05 


— 3.15 ± 0.04 


— 3.25 ± 0.03 


— 3.29 ± 0.03 






Rk 


5.37 ± 0.05 


A A f 1 /~\ O 

4.45 ± 0.02 


= 3.3 


2.60 ± 0.02 


1.95 ± 0.03 


0.95 ± 0.04 


0.51 ± 0.05 


0.32 ± 0.05 






Pk 


— 1.86 ± 0.04 


— 1.23 ± 0.03 


= —1.00 


— 0.89 ± 0.02 


— 0.87 ± 0.02 


— 0.70 ± 0.02 


— 0.61 ± 0.03 


— 0.58 ± 0.03 


3-15 


7270 


v k 


14.36 ± 0.03 


13.84 ± 0.02 


13.02 ± 0.02 


12.74 ± 0.02 


12.36 ± 0.02 


11.88 ± 0.02 


11.53 ± 0.02 


11.46 ± 0.02 






a 'k 


-1.86 ± 0.12 


-2.28 ± 0.09 


-2.66 ± 0.07 


-2.81 ± 0.06 


-2.97 ± 0.05 


-3.15 ± 0.04 


-3.26 ± 0.04 


-3.30 ± 0.04 






Rk 


5.23 ± 0.06 


4.33 ± 0.03 


= 3.3 


2.62 ± 0.02 


2.00 ± 0.03 


0.89 ± 0.05 


0.51 ± 0.06 


0.32 ± 0.06 






Pk 


-1.86 ± 0.05 


-1.25 ± 0.03 


= -1.00 


-0.87 ± 0.03 


-0.86 ± 0.02 


-0.68 ± 0.03 


-0.60 ± 0.03 


0.57 ± 0.03 






~t'k 


0.34 ± 0.22 


0.20 ± 0.18 


-0.08 ± 0.14 


-0.18 ± 0.14 


-0.21 ± 0.14 


-0.26 ± 0.13 


-0.34 ± 0.14 


-0.40 ± 0.13 


Priors 




Rk 


5.05 ± 0.10 


4.31 ± 0.10 


= 3.3 


2.73 ± 0.10 


2.07± 0.10 


0.95 ± 0.10 


0.64 ± 0.10 


0.39 ± 0.10 



(3 k -1.84 ± 0.05 -1.40 ±0.05 =-1.00 -0.84 ± 0.05 -0.70 ± 0.05 -0.52 ± 0.05 -0.38 ± 0.05 0.36 ± 0.05 



Note. — All correlations assume an LMC distance modulus of 18.5 mag. The mean LMC extinction is E(B — V) = 0.10 for MFl, MF2, Model 0, and Model 1, and it is 
E(B - V) = 0.15 for Models 2-15 and 3-15. MFl and MF2 are the Madore & Freedman (1991) PLC relations for the LMC after shifting the period origin to logP = 1-4 from 
1.0. The uncertainties in MFl and MF2 the zero-points may be exaggerated by the shift in the period origin because we did not possess the correlations in the a'^ errors. 
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Table 3. Cepheid Distances & Extinctions 



Galaxy 


Var 


Published 


Model 


Model 1 


Model 2-15 


Model 3-15 


vjaiaxy 


Rn 


L'-t.OA ZC U.1J 


14 37 -U n 17 
14.0/ zc u. ± t 


14 19-1-0 11 


1 A 49 -I- 11 


1 4 ^1 4- fl 1 9 
11.00 in u.j-^ 




Pint km <i _1 ~l 

W(J I Kill b 1 


254 ± 21 


244 ± 18 


240 ± 13 


241 ± 14 


241 ± 14 


LMC 


^ — f^LMC 




= 00 


— o.o 


= 00 


= 0.0 




/ Ft B — V")\ 


0.08 


— o.io 


— 0.10 


— 0.15 


— 0.15 






0.03 




u.uyo zc u.uuu 


07Q -1- 007 
u.u ( y u.uu i 


n f|7Q -1- fl 007 
u.u / y u.uu/ 


SMC 


M — fJ'LM C 


0.4 ± 0.1 


47 4- 02 

U.rt I _1_ U.Uii 


n « + o 112 

U.OO _l_ U . L/Zj 


(1 48 + (1 02 

U.rtO _1_ U.Ui 


*M + o 05 

U.OO _l_ u.uo 




(FAR — V\\ 


0.04 


n 052 + n m n 


037 + 003 

U.UOl _l_ u.uuo 


1 00 + 003 


n 1 + 01 4 

U. -LO 1 _l_ U.U-Lt: 




& E 


0.02 




089 4- 004 

U.UOZ ZC U.UU1 


Ofifi 4- 00^ 
u.uuu zc u.uuo 


ORR 4- OO 1 ^ 
u.uuu zc u.uuo 


M 33 


V — PLM C 


i a _|_ n no 
u. ±1 zc u . uy 


r 09 4- o i ^ 

U.UZr ZC U.-LO 


« 00 4- OR 
u.uu zc u.uu 


f. oo 4- 0^ 
u.uu zc u.uo 


3cc 4. OQ 
u.oo zc u.uy 




I Ft B — 1/1 \ 


n 1 n 4- OQ 
u. ±u zc u . uy 


n ifil + n 0^0 


1 fi9 4- 090 

U. lUi zc u.u^u 


90Q 4- 01 R 
u.Ziuy ziz u.U-Lu 


1 07 4- 030 

U.LUI ZC U.UOU 




G E 






O^ 1 ^ 4- 01 ^ 
u.uoo zn u.u±o 


019 4-0 0^0 

U.U-LZr ZC U.UOU 


01 4- 0^0 

U.U-LU ZC U.UOU 


M 31 


p — f^LMC 


C CM _|_ fl 1 A 


O.OO ZC U. 1U 


^ 77 4- 04 

i). I t ZC U.U4: 


^ 80 4- 03 
o.ou zc u.uo 


^ Q7 + nn c i 
o. y i zc u.uo 




(Ft B — V\\ 


n 10 4-0 13 
u. ± y zc u.io 


1 Q4 4- n 039 


911 4- o 01 9 

U.il 1 Z1Z U.UJ-ii 


9AQ 4-0 010 
u.Zr^iy zc u.uxu 


900 4- 01 ^ 

U.^UU ZC U.U-LO 




& E 






1 Ofi 4- 01 9 

U.1UU ZC U.UJ-Zr 


1 1 4- 01 R 

U..L-LU ZC U.U±U 


1 O^t 4- 01 Q 
U.-LUU zc u.u±y 


IN V_i V_, OUU 


M — fJ-LMC 


0. 1U ZC U.1U 


oil _|_ n 1 4 

O. 1 -L ZC U. 14: 


Q 1 Q 4. 03 
O.IO X u.uo 


8 09 4- 03 
O.UZr zc u.uo 


dropped 




/ Ft B — V")\ 


_n 07 -i- n en 

— u.u/ zc u.uo 


—0 037 4- 047 

— U.UO 1 ZC U.U1 1 


—0 048 4- 01 

U.UIO ZC U.U-LU 


034 4- 008 
U.UOl zc u.uuo 












0^2 + 00Q 

u.uji _l u.uuy 


+ 01 fi 




M 81 


^ — PLMC 


Q Qfl _|_ o 90 
y.ou zc u.ZrU 


q 99 4-0 91 


q 99 4- o 0Q 
y.^Zi zc u.uy 


q 0^ 4- OQ 
y.uo zc u.uy 


q 1Q4-0 11 

y . -LO ZC U. -L _L 




/ Ft B — V")\ 


03 4- 0^ 
u.uo zc u.uo 


07^ 4- 07R 
u.u i o zc u.u t u 


071 4- 097 
U.U / J. ZC u.u^/ 


1 80 4- 037 

U.-LOU ZC U.UOl 


1 83 4- 03R 

U..LOO ZC U.UOU 




°" E 






OQ^i 4- 01 9 
u.uyo zc u.u±Zr 


079 4-0 018 

U.U 1 ^ ZC U.U-LO 


OR 1 ^ 4- 01 7 

U.UUO ZC U.U.L 1 


M 101 


p — PLMC 


1 84 4- 17 

-LU.O^i ZC U.A I 


i o £9 4- 91 

J.U.OZ1 ziz u.^x 


i o 84 4- 1 1 

1U.04 ZIZ U.J--L 


i RR 4- OQ 
lu.uu zc u.uy 


in ^R 4- 11 

LU.UU ZC U.-LJ. 




1 Ft B — 


0.03±?? 


01 ^ 4- 07Q 
u.u±o in u.u f y 


00Q 4- 094 
u.uuy zc u.u^i 


1934-0 033 

U.-L-^O ZC U.UOO 


1 RR 4- 044 

U..LUU ZC U.U'il 










OQ^i 4- 01 
u.uyo zc u.uj-U 


07Q 4-0 01^ 
u.u i y zc u.U-LO 


07Q 4-0 017 
u.u I y zc u.ul / 


IC 4182 


M - Plmc 


9.86 ± 0.09 


9.98 ± 0.23 


10.05 ± 0.08 


9.81 ± 0.11 


9.57 ± 0.14 




(E(B-V)) 


-0.08 ±0.07 


-0.041 ±0.076 


-0.064 ± 0.031 


0.065 ±0.039 


0.153 ±0.056 




<?E 






0.081 ± 0.013 


0.044 ± 0.022 


0.044 ± 0.023 


NGC 5253 


P- - PLMC 


9.58 ±0.10 


9.20 ±0.32 


9.23 ±0.16 


9.07 ±0.13 


8.90 ± 0.15 




(E(B-V)) 


0.02 ±0.29 


0.071 ±0.125 


0.058 ± 0.061 


0.165 ±0.046 


0.209 ±0.051 




<JE 






0.080 ± 0.023 


0.044 ± 0.038 


0.054 ±0.032 


NGC 925 


M - V-LMC 


11.34 ±0.16 


11.33 ±0.18 


11.31 ±0.08 


11.18 ±0.10 


11.23 ± 0.10 




(E(B-V)) 


0.13 ±0.08 


0.136 ±0.067 


0.137 ±0.035 


0.235 ±0.031 


0.244 ±0.034 




VE 






0.102 ± 0.007 


0.069 ±0.010 


0.069 ±0.011 


M 96 


P- - PLMC 


11.82 ±0.16 


11.78 ±0.28 


11.83 ±0.10 


11.64 ±0.10 


11.92 ± 0.15 




(E(B-V)) 


0.06 ±0.03 


0.061 ±0.117 


0.043 ± 0.029 


0.159 ±0.043 


0.127 ±0.038 




(T E 






0.093 ± 0.023 


0.076 ± 0.029 


0.067 ±0.030 
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Table 3 — Continued 



Galaxy 


Var 


Published 


Model 


Model 1 


Model 2-15 


Model 3-15 


NGC 3351 


H — P-LMC 


11.51 ±0.19 


11.42 ±0.19 


11.41 ±0.06 


11.27 ±0.08 


11.59 ± 0.14 




(E(B-V)) 


0.12 ±0.02 


0.122 ±0.073 


0.122 ± 0.047 


0.222 ±0.033 


0.176 ±0.032 




&E 






0.122 ± 0.013 


0.105 ±0.014 


0.104 ±0.014 


NGC 4536 


M — MiMC 


12.60 ±0.13 


12.47 ±0.22 


12.35 ±0.11 


12.22 ±0.09 


12.11 ± 0.12 




(E(B-V)) 


0.04 ± 0.04 


0.077 ±0.079 


0.117 ±0.032 


0.215 ±0.035 


0.261 ±0.042 










0.083 ± 0.011 


0.065 ±0.015 


0.064 ±0.015 


M 100 


M — P-LMC 


12.54 ±0.17 


12.49 ±0.20 


12.50 ±0.10 


12.32 ±0.09 


12.60 ±0.15 




(E(B-V)) 


0.10 ±0.06 


0.075 ± 0.072 


0.071 ± 0.020 


0.182 ±0.031 


0.143 ±0.029 










0.117 ±0.011 


0.098 ±0.012 


0.097 ±0.012 


NGC 4496A 


M - /»LMC 


12.53 ±0.14 


12.41 ±0.20 


12.43 ±0.11 


12.27 ±0.10 


12.17 ±0.13 




(E(B-V)) 


0.03 ± 0.04 


0.061 ±0.074 


0.056 ± 0.028 


0.162 ±0.037 


0.207 ± 0.044 




°E 






0.094 ± 0.008 


0.066 ±0.011 


0.066 ±0.011 


NGC 4639 


M — V-LMC 


13.53 ±0.22 


13.61 ±0.32 


13.59 ±0.21 


13.24 ±0.11 


13.28 ± 0.13 




(E(B-V)) 


??±?? 


-0.083 ±0.114 


-0.078 ± 0.072 


0.091 ±0.035 


0.093 ± 0.039 










0.112 ± 0.021 


0.097 ±0.030 


0.098 ± 0.031 



Note. — The published extinction estimates for M 33, M 31, IC 4182, NGC 5253, NGC 4536, NGC 4496A, and NGC 4639 are 
not directly comparable to our estimates because the authors used different extinction models. A ? indicates that no estimate 
was included in the published results. The value for the solar radius Rq is the consensus value from Reid (1993), and the value 
for the circular velocity of the sun 0o combines the Ro estimate with the proper motion of Sgr A* (Backer 1996). The relative 
distance of the LMC and the SMC and the extinction values are from the review by Westerlund (1990), although Jacoby et 
al. (1990) find higher values of 0.13 for the LMC and 0.06 for the SMC based on Balmer decrements in planetary nebulae. 
Our distances and the Westerlund (1990) value for the SMC are larger than the Caldwell & Laney (1991) values because of 
differences in defining the Cloud centers. The distance errors are very strongly correlated and cannot be treated as 
independent, random uncertainties! 
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Table 4. Comparison Distances 



Galaxy 


EPM 


SNIa 


SBF 


SNIa/MLCS 








p /^M31 


M — MS-/V1981S 


M 31 






= 0±0.06 




M 81 






3.26 ±0.08 




M 101 










IC 4182 




q 77 +0.44 

y -' '-0.55 






NGC 5253 


(9.56 ±0.2) 


9.5llS:g 


3.65 ±0.10 


-3.17 ±0.09 


NGC 925 


(11.63±°;lt) 




5.53 ±0.09 




M 96 


(12.38t?; 9 3 ?) 




(5.76 ±0.06) 




NGC 3351 


(12.38t?; 9 3 ?) 




(5.76 ± 0.06) 




NGC 4536 




12.89l° ;« 




= 0±0.07 


M 100 


12-3718$ 








NGC 4639 


i3.ooig:£ 






0.66 ±0.14 



Note. — Absolute distances for the Cepheid galaxies based on the expanding photosphere method (EPM, 
Eastman, Schmidt & Kirshncr 1996) or physical models of Type la supernovae (SNIa, Hofiich & Khokhlov 
1996), and relative distances based on the surface brightness fluctuation method (SBF, Tonry et al. 1996) 
and the MLCS method for Type la supernovae (SNIa/MLCS, Riess et al. 1996). The SBF distances are 
relative to M31, and the MLCS distances are relative to SN 1981B in NGC 4536. Values without parenthesis 
are distances to the Cepheid galaxy, while values in parenthesis are for galaxies in the same group. 



